Variational approach to the excitonic phase transition in graphene 
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We analyze the Coulomb interacting problem in undoped graphene layers by using an excitonic 
variational ansatz. By minimizing the energy, we derive a gap equation which reproduces and 
extends known results. We show that a full treatment of the exchange term, which includes the 
renormalization of the Fermi velocity, tends to suppress the phase transition by increasing the critical 
coupling at which the excitonic instability takes place. 
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Introduction. The role of Coulomb interactions in 
the low-energy regime of undoped graphene layers has 
arisen great interest and still remains somewhat contro- 
versial. This is rooted to the poor screening properties 
of graphene, a peculiarity that can be traced back to the 
linear density of states of the low-energy theory, which 
vanishes at the Dirac point. An early weak-coupling 
analysis of the problem, based on the Renormalization 
Group (RG) method, showed that the Coulomb interac- 
tion is marginally irrelevant, flowing to a non-interacting 
fixed point i This picture can be rigorously justified at 
all couplings in the limit of a large number of electron 
flavors^ In that scenario, undoped graphene layers would 
behave mostly as a non-interacting system of electrons, 
with minor traces of interactions reflected in the lifetime 
of quasiparticles 3 - and in a logarithmic renormalization 
of the Fermi velocity. 1 Remarkably, this picture seems to 
match reasonably well with current experimental data4^— 

The relative strength of the Coulomb interaction mea- 
sured, as compared with the kinetic energy of the elec- 
trons, is ruled by the dimensionless coupling constant 
9o = e 2 /ei>i?, where e is the electron charge, vf the Fermi 
velocity, and e the dielectric constant of the medium in 
which graphene is embedded. We use the subscript to 
denote unscreened values of the coupling, see below. For 
samples in vacuum, e 2 ~ 14.4 eV A, and vp ~ 6.3 eV 
A, so that go ~ 2.3. Density Functional Theory 7 ^ give a 
value for the screened coupling in the range g ~ 0.5 — 2. 
This puts graphene in the intermediate coupling regime, 
and hence the validity of the weak-coupling analysis re- 
lies on the absence of a strong coupling fixed point in 
the RG transformation. Indications of such a fixed point 
have been found by extending the weak-coupling RG to 
higher orders in the coupling constant expansion. 9 The 
experimental data would still be compatible with this 
strong-coupling scenario if current setups had graphene 
sufficiently isolated from the environment and could op- 
erate with perfectly neutral samples. 

The possibility of phases beyond the reach of pertur- 
bative or weak coupling renormalization group methods 
in undoped graphene has been explored in the literature 
by using different approaches. The main candidate for 
a strong coupling phase is an excitonic condensate, in 



which electron and holes bind together opening up a 
gap in the density of states and rendering the system 
insulating. The mechanism responsible for this phase 
would be the gain in exchange energy arising from the 
long range Coulomb interaction. A gap equation for this 
transition has been derived within the Dyson- Schwinger 
formalism^ and different solutions of this equation yield 
(unscreened) critical couplings for the phase transition 
around an unscreened coupling g$ c ~ 1 — 2.— ~— This 
scheme is equivalent to the summation of a class of dia- 
grams, and can be considered an extension of weak cou- 
pling approaches. Montecarlo calculations in the lattice 
have been carried out to analyze this problem , 13 ' 14 find- 
ing an insulating phase above go c ~ 1.11 and go c ~ 1.66, 
respectively. A phase transition beyond a certain cou- 
pling can also be found for short range interactions in 
the half filled honeycomb lattice, 15 although the critical 
coupling takes the model beyond the regime where the 
approximation of the electronic bands by the Dirac equa- 
tion is valid. Finally, the study of the two-body problem 
in graphene, with Coulomb interactions, leads to a re- 
markable instability of the wave-function for a critical 
coupling go c — l, 16 that might underlie the eventual for- 
mation of excitons. 

In this paper we investigate the possibility of an ex- 
citonic strong coupling phase in undoped graphene by 
using a variational ansatz. The method used here can 
be extended in a straightforward way to finite temper- 
atures or to finite carrier concentrations. We derive a 
gap equation similar to that obtained in the previous 
literature ; 10 ' 12 but with the inclusion of a the renormal- 
ization of the Fermi velocity. By analyzing numerically 
and analytically the resulting gap equation, we find that 
the latter produces a supprcsion of the phase transition 
by increasing the critical coupling. 

The model and variational ansatz. As we have men- 
tioned, we will address the problem of an undoped 
graphene sample with Coulomb interactions in the low- 
energy regime, where the electron motion is described by 
the Dirac equation. The Hamiltonian for this problem 
reads: 
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FIG. 1: Diagrams included in the ground state energy within 
the variational ansatz. a) Hartree term, which is zero, b) 
Exchange term, which is the dominant one. 



where s = ± refers to the upper and lower cones, re- 
spectively. In our model, the spin and valley degrees of 
freedom are considered only as extra degeneracies in the 
number of fermions. The Coulomb potential is given by 
V q = 2ng J ' q, where g is the dimensionless coupling con- 
stant introduced above. We assume that g includes con- 
tributions from static screening^S. In the RPA approx- 
imation, transitions between the valence and conduction 
bands lead to a momentum independent dielectric con- 
stant, which can be incorporated in a straightforward 
into this formalism. Alternatively, one can view this ap- 
proach as the leading approximation in the limit N — > oo, 
where N is the number of fermion flavors 17 . 

Our goal is to analyze the ground state of this Hamil- 
tonian by using a variational ansatz which includes the 
possibility of pairing between electron and holes. Such 
an ansatz was proposed to study excitons formation in 
semiconductors, 18 and is reminiscent of the ansatz used 
in the BCS theory of superconductivity: 
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The ansatz contains a coherent superposition of states 
with a different number of electron- hole pairs. Here \D) 
stands for the filled Dirac sea, and u k and v k are vari- 
ational parameters to be determined by minimizing the 
ground state energy. Without loss of generality, they are 
taken real. Notice that they are not independent, since 
the normalization of the wave function imposes the con- 
straint: 
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Derivation of the gap equation. Following the lines of a 
typical variational calculation, the energy of the ansatz is 
evaluated by projecting the Hamiltonian into this state. 
It has two contributions, the Hartree and the exchange 
one, as shown in Fig. [TJ The Hartree contribution is 
zero, by virtue of the normal ordering of the Hamiltonian 
with respect to the Dirac sea, which is physically related 
to the neutrality of charge of the global system. The 
dominant contribution comes from the exchange energy, 
which includes terms with a momentum transfer of q = 



k' — k. The projected Hamiltonian reads then: 
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where we have used the normal-ordered Hamiltonian in 
order to carry out the calculation. The extreme condition 
must be imposed respecting the normalization constraint. 
The result gives the following equation: 
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This equation can be simplified by introducing the fol- 
lowing parameters: 
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The first equation is the self-energy insertion to the elec- 
tron propagator, which adds to the linear term coming 
from the non-interacting dispersion relation and repre- 
sents a renormalization of the Fermi velocity. The second 
equation introduces A^, which can be identified with the 
gap that arises in the electronic spectrum when excitons 
are formed. This is clearly expressed in the third equa- 
tion, which gives the dispersion relation of Bogoliubov 
quasiparticles in the excitonic condensante. 

In terms of these new parameters, the solution to the 
variational problem reads: 
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By plugging these expressions into the equation for the 
gap, Eq. ([7]), we get a self-consistent integral equation, 
namely: 
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As we have already mentioned, a similar gap equation 
has been already found by using the Schwinger-Dyson 
formalism. 10 

Further insight can be obtained by carrying out the 
angular integral while keeping the lowest order terms in 
a Legendre polynomial expansion of the Coulomb inter- 
action Vjk'_ kl • This yields a simplified integral equation 
in the continuum limit of the problem: 
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where we have introduced the following kernel: 
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The main feature of this gap equation, as compared with 
previous approaches, is the inclusion of the exchange cor- 
rection to the free electron dispersion relation, Eq. (J6j) . 

Analysis of the gap equation. In order to extract in- 
formation from the gap equation, we make the follow- 
ing assumption! 10 i 19 the dominant contribution to the 
gap equation corresponds to the region fc 3> A*, where 
A* = Aa* (with vp = 1). This allows us to make 
E q ~ £q and write A* as the lower limit of the integral: 
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By using the same type of reasoning, an expression for 
the £& can be derived, which only retains the leading, 
most divergent terms (and valid for £& > A*): 
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As mentioned above, this is actually the renormalization 
of the Fermi velocity that arises from a RG analysis, 1 
which has been so far neglected in the literature on the 
excitonic condensation. We will see shortly that this log- 
arithmic correction plays a crucial role in the analysis of 
the gap equation. 

Let us transform the integral equation into a differen- 
tial equation, namely: 
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which has the form of a 3D radial Schrodingcr equation 
with a potential g(k). The latter is the running cou- 
pling constant in the RG sense, which has appeared in 
a natural way from the exchange correction to the linear 
dispersion relation. It reads: 
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The differential equation (I17p must be supplemented with 
boundary conditions that are also derived from Eq. (|T5|) : 

fc 2 A' fe | fe=A . = 
(fcA' fe + A fc )| fc =A = 



(19) 

(20) 



The first one is the infrared condition, since it is eval- 
uated at the gap A* -c A, while the second one is the 
ultraviolet one, evaluated at the cutoff. 

Adiabatic solution. A preliminary study of Eqs. (|17[) - 
([20]) can be made by assuming that g(k) varies slowly 
enough for an adiabatic approximation to be reasonable. 
Noting that the case of a constant potential g(k) = g 
admits an exact solution of the form 
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the adiabatic solution can be found to be, following Ref. 
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where f(k) = J" %-\l y{k') — \ ■ Implementation of 
the boundary conditions (TTT)l) - (|20p yields the quantiza- 
tion rule 
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where n is a positive integer and 5k = arctan y/4g(k) — 1. 
The goal is to solve for A* , a nonzero value meaning that 
there is an excitonic instability. The condition g(k) > 1/4 
for all values of k, leads to the requirement A* > A min — 
Ae~ 8 ( 1 ~ 1 / 4g \ We find that a nonzero, real solution of 
Eq. (f2"2")l (with n = 1) satisfying A* > A min , exists for 
g greater than a critical value g c ~ 0.5, which marks the 
onset of the excitonic instability. 

Numerical solution. We further check the previous 
analysis by numerically solving Eq. (|17[) with the bound- 
ary conditions (fi"9|) - ([20|) . The results are shown in Fig. [2] 
We find solutions for g > <? c ~ 0.59, in reasonable agree- 
ment with the adiabatic approximation. The asymptotic 
limit Afc ~ 1/vfc [see Eq. (J21I) ] is only clearly visible 
for g ~ g c and fc « A. A detailed analysis of the region 
where (g — g c ) — > + suggests that A* ex (g — g c ), see 
Fig. [3] For comparison, we also show the numerical re- 
sults obtained by neglecting the renormalization of the 
Fermi velocity. They reproduce correctly the main fea- 
tures found in analytical studies, namely, g c = 1/4 and 
A* ~ e - A /Va-9c^ w here A is a constant^— 

Conclusions. We have analyzed the problem of 
Coulomb interactions in undoped graphene by using a 
variational ansatz that includes the possibility of exci- 
ton formation. Our approach can be readily extended 
to other two- and three-dimensional materials, as well as 
to finite temperatures and carrier concentrations. 20 It al- 
lows us to calculate the total free energy, which can be 
compared to that of other broken symmetry phases. 

Our variational analysis reproduces the main features 
of the excitonic transition in graphene. 10 In addition, we 
find that a renormalization of the Fermi velocity is a nat- 
ural by-product of the variational treatment. The result- 
ing change in one particle energies leads to a cancela- 
tion of the leading divergences with trigger the excitonic 
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FIG. 2: (Color online). Left: Dependence of A* on g ob- 
tained by numerically integrating Eq. (|17p with the boundary 
conditions in Eqs. (|19[) -(120 [ ). Right: Dependence of Ak on 
k for g — 0.6,0.65,0.7,0.75 and 0.8 (from top to bottom). 
The thick black line shows the position of the maxima of A^ , 
which give the value of A*, see Eq. (|19|l . 
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FIG. 3: (Color online). Details of the dependence of A* on 
g near the transition. The inset on the right graph shows 
the results from numerically solving Eq. ()17|) without velocity 
renormalization [i.e. assuming g(k) = g in Eq. (I18|l ]. 



transition. A similar effect is observed in the analysis of 
the excitonic transition due to short range interactions in 
graphene bilayers . 21 ' 22 Our variational analysis leads to 
a critical coupling g c ~ 0.59, which is about a factor two 
larger than the critical coupling obtained neglecting the 
Fermi velocity renormalization, g c = l/4 i 10 i 12 i 23 i 24 The 
renormalization of g to two loops gives a transition at 
g c ~ 0.83. 9 



In suspended graphene, we can assume that the effec- 
tive value of g is modified solely by internal screening^iiS 
so that g = go/(l + Nirgo/8), where N — 4 is the 
number of electron flavors, and go — e 2 /vp is the 
bare coupling constant. Then, an upper bound of g is 
ffmax = \im e 2/ VF _ ¥OCl g(e 2 /v F ) = 2/n ~ 0.64, which lies 
slightly above the value of g c obtained with our varia- 
tional ansatz. For the realistic value of go = 2.3 we ob- 
tain g = 0.50, which is below the critical value g c ~ 0.59 
which we have found by solving Eq. (JTTJ) numerically. 25 



Acknowledgements. We acknowledge financial sup- 
port by MICINN (Spain) through grants FIS2007-65723, 
FIS2008-00124and CONSOLIDER CSD2007-00010, and 
by the Comunidad de Madrid, through NANOBIOMAG 
and MICROSERES. We also acknowledge useful discus- 
sions with J. Gonzalez, M. I. Katsnelson, M. Polini, and 
I. Zapata. 



1 J. Gonzalez, F. Guinea, and M. A. H. Vozmediano, Nucl. 
Phys. B 424, 595 (1994). 

2 J. Gonzalez, F. Guinea, and M. A. H. Vozmediano, Phys. 
Rev. B 59, R2474 (1999). 

3 J. Gonzalez, F. Guinea, and M. A. H. Vozmediano, Phys. 
Rev. Lett. 77, 3589 (1996). 

4 Z. Q. Li, E. A. Henriksen, Z. Jiang, Z. Hao, M. C. Martin, 
P. Kim, H. L. Stormer, and D. N. Basov, Nature Physics 
4, 532 (2008). 

5 R. Nair, P. Blake, A. Grigorenko, K. Novoselov, T. Booth, 
T. Stauber, N. Peres, and A. Geim, Science 320, 1308 
(2008). 

6 G. Li, A. Luican, and E. Y. Andrei, Physical Review Let- 
ters 102, 176804 (2009). 

7 M. Polini, A. Tomadin, R. Asgari, and A. MacDonald, 
Phys. Rev. B 78, 115426 (2008). 

8 M. van Schilfgaarde and M. I. Katsnelson (2010), 
arXiv: 1006.2426. 

9 O. Vafek and M. J. Case, Phys. Rev. B 77, 033410 (2008). 

D. V. Khveshchenko, Phys. Rev. Lett. 87, 246802 (2001). 

1 D. Khveshchenko, J. Phys.:Condens. Matter 21, 075303 
(2009). 

2 O. V. Gamayun, E. V. Gorbar, and V. P. Gusynin, Phys. 
Rev. B 81, 075429 (2010). 

3 J. E. Drut and T. A. Lahde, Phys. Rev. Lett. 102, 026802 
(2009). 



14 W. Armour, S. Hands, and C. Strouthos, Phys. Rev. B 81, 
125105 (2010). 

15 S. Sorella and E. Tosatti, Europhys. Lett. 19, 699 (1992). 

16 J.Sabio, F. Sols, and F. Guinea, Phys. Rev. B 81, 045428 
(2009). 

17 M. S. Foster and I. L. Aleiner, Phys. Rev. B 77, 195413 
(2008). 

18 C. Comte and P. Nozieres, J. Physique 43, 1069 (1982). 

19 E. Gorbar, V. Gusynin, V. Miransky, and I. Shovkovy, 
Phys. Rev. B 66, 045108 (2002). 

20 P. Hosur, S. Ryu, and A. Vishwanath, Phys. Rev. B 81, 
045120 (2010).' 

21 R. Nandkishore and L. Levitov (2010), arXiv:1006.1136. 

22 Y. Lemonik, I. L. Aleiner, C. Toke, and V. I. Fal'ko (2010), 
arXiv: 1006. 1399. 

23 E. V. Gorbar, V. P. Gusynin, V. A. Miransky, and I. A. 
Shovkovy, Phys. Rev. B 66, 045108 (2002). 

24 D. V. Khveshchenko (2008), arXiv:0807.0676. 

25 Montecarlo results suggest a critical coupling in the range 
of expected values in suspended samples 13 , which is smaller 
than the value predicted by our variational ansatz. The 
origin of this discrepancy is unclear. It could arise from 
differences in the regularization of the short range part of 
the interaction, or in the treatment of screening at high 
frequencies. 



